Schematic models for dynamic yielding of sheared colloidal glasses 

Matthias Fuchs* and Michael E. Gates 
Department of Physics and Astronomy, The University of Edinburgh, 
JCMB King's Buildings, Edinburgh EH9 3JZ, GB 
(Dated: February 1, 2008) 

Abstract 

The nonlinear rheological properties of dense suspensions are discussed within simplified mod- 
els, suggested by a recent first principles approach to the model of Brownian particles in a 
constant-velocity-gradient solvent flow. Shear thinning of colloidal fluids and dynamical yielding 
of colloidal glasses arise from a competition between a slowing down of structural relaxation, 
because of particle interactions, and enhanced decorrelation of fluctuations, caused by the shear 
advection of density fluctuations. A mode coupling approach is developed to explore the shear- 
induced suppression of particle caging and the resulting speed-up of the structural relaxation. 
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I. INTRODUCTION 



Soft materials, such as particle dispersions, exhibit a wide range of rheological properties. 
While dilute colloids flow with a viscosity only slightly higher than that of the solvent, 
concentrated dispersions behave as weak amorphous solids. For intermediate concentrations, 
one generally observes, upon increasing the external shear rate, first, a strong decrease of 
the dispersion viscosity ("shear thinning"), and then an (often dramatic) increase of the 
viscosity ("shear thickening" )0'i. 

While shear-induced crystallization of particle suspensions causes a marked decrease of 
the viscosity, shear thinning is not always accompanied by a flow- induced orderingii'ii, 
but appears connected more generally to a decrease of the Brownian contribution to the 
stressi'i. In concentrated suspensions of polydisperse colloidal particles the structure re- 
mains amorphous during the application of shear but still exhibits shear thinning or yield 
behaviouii0. 

Detailed light scattering studies of quiescent colloidal hard sphere suspensions0'0'E10 
have identified a slowing down of the structural relaxation as the origin of the solid-like 
behavior at high concentrations. Operationally, a transition to an amorphous solid or glass 
can be defined when the structural relaxation time increases beyond the experimental ob- 
servation time, and various concomittant signatures of metastability have been observed. 
The resulting amorphous solids, even though their life-time may be limited by aging and 
crystallisation processes, nevertheless were found to possess a well defined average arrested 
structure. On the theoretical side, predictions for the glassy structure and the slow-down of 
the structural relaxation have been obtained within the mode coupling theory (MOT), which 
describes an idealized glass transition scenario with a divergent structural relaxation time 
at a critical concentration (or temperature )lli@0. Comparisons of theory and experimental 
data have shown agreement on the 30% relative error level§@'i'0'§. 

Considering that arrested systems like coUidal glasses often are redispersed by shaking 
or stirring the sample, the influence of external shear strain or stress on glasses obviously 
is of interest. Because glass formation is connected to a growing internal relaxation time, 
an important aspect of the imposition of external driving is the introduction of new time 
scales. In the case of steady shearing, this aspect is generally discussed as interplay of shear- 
induced and Brownian motion. Because, without shear, the glassy system fails to reach 
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equilibrium for long times, it is unclear to which stationary non-equilibrium state shear 
motion leads, and how this non-equilibrium state approaches the equilibrium one in the 
limit of vanishing shear rate. Various phenomenological models ( "constitutive equations" ) 
either describe a yield-stress discontinuity (Bingham plastic, Hershel-Bulkeley law), or a 
power-law approach (power-law fluid)0, while recent glass theories predict the formeJil, the 
lattei@, or a transition from one to the other with concentration or temperatureEiJ. 

In order to gain more insight into the yielding of colloidal glasses and the nonlinear rhe- 
ology of colloidal fluids, we extend here the analysis of the microscopic approach which we 
recently presented^. We derive from it simplified models, aiming to bring out the quali- 
tative and universal features. A simple case of non-linear rheology is studied: a system of 
Brownian particles in a prescribed steady shear solvent flow with constant velocity gradient. 
While hydrodynamic interactions and fluctuations in the velocity profile are thus neglected 
from the outset, this (microscopic) model has the advantage that the equation of motion 
for the temporal evolution of the many-particle distribution function can be written down 
exactlyi'ii. Our theoretical development can therefore be crucially tested by Brownian dy- 
namics simulations, like Ref.i, and constitutes a first microscopic approach to real glassy 
colloidal suspensions. The properties of this microscopic model have been worked out for 
low densitieJ^lH and it provides the starting point for various (approximate) theories for 
intermediate concentrations^. 



II. MICROSCOPIC APPROACH 



The effect of a constant uniform shear rate 7 on the particle dynamics is measured by the 
Peclet number!, Pcq = 'jd'^/DQ, formed with the bare diffusion coefficient Dq of a particle of 
diameter d. If the quiescent systems exhibits a much longer time scale r, as do dense colloidal 
suspensions where r is the final or structural relaxation time, then a second, "dressed" Peclet 
(or Weissenberg) number, Pe = 7r, can be defined. This characterizes the influence of shear 
on the structural relaxation and increases without bound at the glass transition, even while 
Peo < 1. InRef.li, we argued that the competition of structural rearrangement and shearing 
that arises when Pe > 1 ^ Peo dominates the non-linear rheology of colloids near the glass 
transition. 

While the time scale ratios Peo and Pe appear quite generally, the physical mechanisms 
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active when Pe > 1 may be quite different for different classes of soft matter. As we argued 
in the introduction, flow- induced ordering shall here be assumed absent. In the quiescent 
dispersion, the structural relaxation is dominated by (potential) particle interactions which 
either cag or bondii a central particle among its neighbours. Either mechanism leads to 
a slowing down of particle rearrangements accompanied by growing memory effects. The 
former process ("cage-effect") is driven by the local order as measured in the height of the 
principal peak of the static structure factor, 5*^^, and leads to a prolonged decay-time of 
especially this density mode; note that its wavevector is inversely related to the average 
particle spacing. The decay time of this dominant cage mode sets the structural relax- 
ation time r, which in the following shall be defined by ^^^(t/r) = 0.1, where $g(t) is the 
normalized intermediate scattering function. 

It is important to realize that during the time interval so defined a single particle has 
diffused (in direction x, say) a fraction of its size onlyHili^: (Ax^(r)) ^ 0.1(P. In the 
sheared system at Pe = 1, therefore, it is not true that kinetic flow of the particle with the 
solvent ("Taylor dispersion" which would give (Ax^(t)) = ^Dq'jH^ along the flow directionll) 
has displaced the particles relative to each other during the time interval t = t. This would 
imply a rapid destruction of the cage already for Pe < 1. Rather, the effect of shear on 
the structural relaxation has to be sought in its effect on the cage mode, i.e., the collective 
density fluctuations with wavevector qp. Because the steady state structure factor differs 
only smoothlyi'Sim from the quiescent Sg around q ^ Qp for Pcq ^ 1, the effect of shear 
cannot lie in a destruction of the steady state local order. This would require larger Peclet 
numbers and contradict the notion that an infinitesimal steady shear rate melts a glass. 
Rather, the effect needs to arise from a shear-induced decorrelation of the memory built up 
in the collective cage mode. The likely mechanism thus appears to be the shear advection 
of fluctuationsiiS 

Figure |l| sketches, neglecting Brownian motion, the advection of a fluctuation with initial 
wavevector (27r/A2;, 0, 0), where x points along the flow direction, into one with wavevector 
{2ti /\x,2t{ l\y{t),Q), where \y{t) = X^/iji) at later time t. Clearly, for any initial fluctu- 
ation with A^. 7^ 0, the wavelength in the gradient direction, Ay(t), will decrease for large 
times. Brownian particle motions (assisted by the interaction forces) "smear out" the fluctu- 
ation with time and cause the decay of the corresponding correlator. Because the advected 
wavelength decreases upon shearing, smaller and smaller motions can cause the fluctuation 
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FIG. 1: Advection by steady shear of a fluctuation in j;-direction with wavelength A^^ at t = 0. At 
later time t, its wavelength in y-direction obeys: Xx/\ = Ax/ Ay = 7t. 

to decay. We presume that the shear advection of density fluctuations with wavevector 
leads to a competition between the non-linear feedback mechanism of the cage effect and 
the shear-induced decorrelation, and that this competition determines the nonlinear rheol- 
ogy of concentrated particle dispersions. Obviously, this competition involves a cooperative 
rearrangement of the (finite number of correlated) particles forming the distorted cage. Any 
microscopic theory like ours therefore requires severe uncontrolled approximations because 
no conceptually well-controlled approximation scheme appears suited to this problem. Thus, 
we build upon MOT (rather than, e.g., approaches starting with uncorrelated binary inter- 
actions, or based on coarse-graining procedures). 



A. Structural dynamics under shear 

We consider a suspension of Brownian particles, with density n = N/V, described by 
the Smoluchowski equation without hydrodynamic interactions. From time t = on, a flow 
is imposed in the solvent, which points in the x- and increases linearly along the y-direction: 
v(r, t) = K,r9(t) = jyexO{t) (where k, is the shear rate tensor, Kij = 'j6ixSjy, and 9(t) the 
step function). Note that we neglect deviations from the imposed linear flow proflle and 
thus cannot capture various shear-banding and other layering phenomena. (The latter may 
or may not arise in real experiments where only the stress, and not the velocity gradient, 
can be assumed constant across a sheared planar sample in steady state.) For this situation, 
the Smoluchowski equation for the particle distribution function is easily formulatedlii. 
But, with shear, its stationary solution is not known in general (except for some results 



at vanishing particle concentrationii'il), and thus steady state quantities and correlation 
functions are not available. 

Recently, Lionberger and Russel have made progress in the semidilute concentration 
regime by transferring liquid state approaches to the steady state situation; see Ref.ii and 
works cited there. Yet, close to the glass transition where the quiescent system develops 
divergent time scales, a liquid state approach does not capture the falling out of equilibrium 
of the system, and thus (presumably) cannot handle the transition to dynamic yielding 
of a metastable solid. In order to capture the inherent long time scales we instead have 
suggested^! to monitor the transient fluctuations of the suspension after turning on the 
solvent flow field at time t = 0. Yet this problem too cannot be solved exactly, and requires 
approximations. First, the relevant variables whose transient dynamics shall be monitored, 
need to be chosen. Then equations of motion for these variables need to be formulated. To 
make progress we build upon the insights into the quiescent system provided by the MCT 
and generalize it to the out-of-equilibrium situation. 

Colloidal suspensions upon densification exhibit a slowing down of the structural relax- 
ation (particle rearrangements) and, to describe it, density fluctuations with wavevector q 
enter as natural variables, = Y^iLi g^^^'- Their magnitude is measured by the equilibrium 
structure factor Sq = {g^gq}/N, which changes rather little, while their dynamics slows down 
dramatically. We follow MCT in considering the set of density fluctuations as a set of slow 
variables. Thus, arbitrary steady-state expectation values are obtained from determining 
the overlap of the relevant quantities with density fluctuations. This requires us to find the 
transient density fluctuations. They shall be determined from a closed set of equations of 
motion for the intermediate scattering functions, which is obtained by performing a mode 
coupling approximation. 

With these approximations, various steady state quantities, like the thermodynamic shear 
stressii cr = (o'xy/V)^'^'^ can be calculated. (Note that steady state averages are abbreviated 
by {. . .)^'^\ while equilibrium ones without shear are given by {■■■)■) The transverse stress 
also provides the shear viscosity which follows as 77(7) = rjoo+cr/'y, where the solvent viscosity 
is denoted rjoo- Our flnal microscopic expression for the steady state stress a is found to be: 

.(7) ^ ^^T'y f klkykyjt) ^fcgfc(t) 2, 



with t the time since switch-on, /c^T the thermal energy, and S*^ = dSk/dk. The tran- 



sient density fluctuations are given by $k(^) = (^'k(t)(^)^'k)/(-^5'fc) and are normalized 
by the quiescent Sk- Because of shear advection, time dependent wavevectors k{t) = 
^Jk'^ + 2kxky'jt + fc^7^t^ appear; a fluctuation with wavevector k at the start of shearing 
t = has a flnite overlap with a fluctuation at the advected wavevector k(t) = k + k/tt at 
the later time t > 0. 

In summary, steady state quantities shall be determined by considering the structural 
relaxation under shearing and "integrating through the transient dynamics" . Because Peo <C 
1, we expect ordering or layering transitions to be absenti; and as hydrodynamic interactions 
are presumed to play a subordinate role during the structural relaxation! we neglect these 
too, focusing solely on the Brownian contribution to the transverse (shear) stress. A major 
approximation within our approach entails the elimination of explicit particle forces in favour 
of the quiescent-state structure factor 5*^ (the only input in our theory and taken to be 
known). This is a near-equilibrium assumption that is formally uncontrolled but motivated, 
at least in part, by the smallness of Pcq. Because we approximate nonlinear couplings under 
shear using equilibrium averages, we require the system to remain "close to equilibrium" in 
some sense. We leave to future work any attempt to make this sense more precise and, if 
possible, to relax the assumption involved. 



B. Glass stability analysis 

Little can be gained in our approach without specifying the dynamics of the transient 
intermediate scattering functions which describe how the equilibrium structure changes with 
time into that of a steadily sheared state. Obviously, any uncontrolled approximation (like 
mode coupling), which we are now forced to perform to obtain the equations of motion for 
the transient fluctuations, can introduce errors of unknown quality and magnitude into our 
results. But while the quantitative accuracy of the equations we propose in Ref.i has not 
been tested yet, there are qualitative conclusions which can be drawn from the structure of 
the equations, which are rather independent of the microscopic details of the approximation. 
Thus, in order to test our basic approach, these more universal aspects are of central interest 
and should be the ones chosen for initial comparison to experiments or simulations. We will 
consider only these aspects here, but in later sections two simple models are used to study 
how far the universal aspects dominate the model-dependent results. 
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The universal results which follow from our approach are connected to the stability equa- 
tion of a (quiescently) arrested glassy structure which gets melted away by the imposed 
shear. From an expansion of the transient density fluctuations around the initially arrested 
structure close to the glass transition we get the bifurcation equations which describe how 
the localization driven by the cage effect competes with the fiuidization induced by shearing. 
The derivation of a generalized "factorization theorem" and the so-called /5-scaling equation, 
proceeds by a straightforward perturbation calculation which is given in Appendix A. Close 
to the bifurcation, the transient density fluctuations are given by one function Q{t), which 
depends on control parameters, and determines the dynamics on all length scales: 

%{t) = f^ + h,g{t). (2) 

The numbers describe the glassy structure at the instability and the critical amplitude 
hq is connected to the cage-breaking particle motion; both retain their definition from the 
unsheared situation. The function Q{t) contains the essential non-linearities of the bifurca- 
tion dynamics which arise from the physical feedback mechanism (the cage effect) and the 
shear disturbance. It depends on a few material parameters only, and follows from 

e - c(^) (7 ty + X g\t) = j^jl dt' g{t - 1') g{t') , (3) 

where £, A and c^^-* are defined in Appendix A. The so-called separation parameter e measures 
the distance to the transition and the exponent parameter A determines power-law exponents 
resulting from Eq. (^, and is known for some systemsEl. The shear-related parameter c*-^^ 
is a (relatively unimportant) number of order unity, which sets the scale for 7 and could be 
absorbed into an effective shear rate y/c^ 7. 

The two derived equations, Eqs. (0,0), describe an expansion around the transition 
point between a non-Newtonian fluid and a yielding solid (within our approach, this ap- 
plies whatever microscopic model is chosen) where divergent relaxation times arise from 
self-consistently calculated memory effects and compete with an externally imposed (shear- 
driven) loss of memory. Corrections of higher order in the small quantities {e,'j,Q) are 
neg lected in Eqs. (|,D; see Ref.ii for a calculation of the leading corrections and for back- 
ground on Eq. (0) at 7 = 0. In the following sections we will discuss two simplified models 
which, on the one hand, allow a more detailed analysis, and on the other hand share the 
universal stability properties derived from Eqs. (@,0). 
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III. MODELS 



We now present two, progressively more simplified, models that provide insights into the 
generic scenario of non-Newtonian flow, shear melting and solid yielding which emerge from 
our approach. 

A. Isotropically sheared hard sphere model (ISHSM) 

On the fully microscopic level of description of a sheared colloidal suspension, kinetic 
flow of the particles with the solvent leads to anisotropic dynamics. Yet, because of the 
strong hindering of the motion at high densities, which leads to a caging of particles, the 
development of the anisotropic "Taylor dispersion" may not yet be important at Pe ~ 1 
so that the motion stays locally isotropic. Recent simulation data of steady state structure 
factors support this consideration and indicate an isotropic distortion of the structure for 
Peo -C 1, while the Weissenberg number Pe is already large lH. A mechanism which operates 
independently from Taylor dispersion arises from the shift of the advected wavevectors with 
time to higher values. As the effective potentials felt by density fluctuations evolve with 
increasing wavevector, this leads to a decrease of friction functions, speed-up of structural 
rearrangements and shear-fiuidization. Therefore, one may hope that a model of isotropically 
sheared hard spheres (ISHSM), which for 7 = exhibits the nonlinear coupling of density 
correlators with wavelength equal to the average particle distance (viz. the "cage-effect"), 
and which, for 7 7^ 0, incorporates shear-advection, may be not too unrealistic. 

Thus, in the ISHSM, the equation of motion for the density fluctuations at time t after 
starting the shear is approximated by the one of the quiescent system: 

$,(t) + r, |$,(t) + Jdt' m,{t - t') $,(t')} = . (4) 

where $g(t) = dt^g(t), and the initial decay rate is Tg = q^D^/Sq with the single particle 
Brownian diffusion coefficient given by the solvent viscosity, Dq = ksT /{"iTirjood). The 
memory function also is taken from the unsheared situation but now shear advection is 
entered into the overlap of the fluctuating forces with density fluctuations at time t: 

^ E nS^(^) ^^(t) -^Iq-kl W , (5) 
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with 

'n?' S S S 

(t) = ■ ^ ^Ht) + q • P [q ■ k Cfc + q • p Cp] (6) 

where p = q — k, and the length of the advected wavevector is approximated by k{t) = 
k{l + {t^f/Sy/^. The effective potentials are given by the direct correlation functions, 
Cq = {1 — l/Sq)/n. (For background on this model without shear, see Refs.i'Sli.) 

For hard spheres, the quiescent Sq, taken from Percus-Yevick theoryi, depends only 
on the packing fraction 0; discretising the wavevector integrations in Eq. as done in 
Ref.ii, we find that the model's glass transition lies at 0c = 0.51591. Thus e (where e ~ 
3.0(0 — 0c)); and 7 are the only two control parameters determining the rheology. The 
exponent parameter becomes A = 0.74 and c^'^^ ~ 3, while Sq has its peak at Qp = 7/d. 

In the same spirit of incorporating advection as the only effect of shearing, the expression 
for the transverse modulus may be simplified to 



The resulting numerical results for the ISHSM are discussed below (section IV). 



B. Schematic model 

The central features of the equations of motion of the ISHSM are that it reproduces the 
stability equation from the microscopic approach, Eq. (^, and that the vertex V^*^^-* con- 
tains the competition of two effects. First, it increases with increasing particle interactions 
("collisions" or "cage effect") which leads to a non-ergodicity transition in the absence of 
shear, and second, it vanishes with time because of shear-induced decorrelation. Both these 
effects can be captured in an even simpler "schematic" model, which moreover can be made 
to obey Eq. (^) also. This schematic F^g^ model considers one normalized correlator $(t), 
with $(t ^ 0) = 1 — Ft . . ., which obeys a generalized relaxation equation: 

$(t) + F |$(t) + Jdt' m{t - t') <i>(t')| = . (8) 

Again, in the absence of the memory kernel m the dynamics is trivial, $(/!:) = exp(— Ft), 
and a low order polynomial ansatz for m suffices to study the generic schematic model. We 
choose 
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Without shear, this model has been studied extensivelyll^o. Increasing particle caging is 
modeled by increasing coupling parameters fi,f2 > 0, and the only effect of shearing is to 
cause a time dependent decay of the friction kernel m. The system loses memory because 
of shearing. The role of the transport coefficient (viscosity) t] is played by the average 
relaxation time obtained from integrating the correlator, and this also is taken to determine 
the stress: 

a = ^r] = ^{T)=^ dt ^{t) . (10) 
Jo 

For the parameters of the model, values studied in the literature shall be taken. While the 
parameter F just sets the time scale, the two interaction parameters are chosen as f 2 = f f = 2 
and Vi = vl + 6vi, where = f 2 ( -y/4/f 2 — 1 ) ~ 0.828. From an analysis similar to Appendix 

A, the critical non-ergodicity parameter is found as /c = 1 ~ 0.293 and the parameters 

in Eq. (^ emerge as follows: 

A = l-/.- 0.707, ,= ^J!lfl±^ = M^^ c(^) = "^^^- + f- ^ 0.586. (11) 

^ Jc ^ Jc ^ Jc 

The major advantage of this schematic model, besides its numerical simplicity, is that 
it encodes the two physical mechanisms at work in two handy parameters. The one, e, 
parametrises the tendency of the undriven system to arrest, while the other, 7/F, measures 
the loss of memory brought about by shearing. 



IV. RESULTS AND DISCUSSIONS 



Both models described above exhibit a non-ergodicity transition which corresponds to an 
ideal fiuid-to-glass transition within MCT. For 7 = 0, upon smooth changes of the input 
equilibrium state parameters, a fluid with $g(t 00) —>■ turns into an amorphous solid, 
$q(t 00) ^ fq > 0. The fg are called glass form factors and describe the arrested 
structure. While transport coefficients of the fluid, like the viscosity, are connected to the 
longest relaxation time r of the $g(t), elastic constants of the solid, like the transverse elastic 
modulus Goo, are given by the fg. (For some in-depth, albeit somewhat older, discussions 
of the idealized MCT see Refs.00, and for a more recent discussion of experimental tests 
Ref.§.) 
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FIG. 2: Correlators <I>(t) of the schematic model for the separation parameters e and shear 
rates j/T as marked in the inset; curves A to F for e = —£i = — 10^^''^^ and 7/r = 10~" with 
n = 0, 2, . . . , 10, while curve G corresponds to an unsheared fluid 7 = at e = — ei (a filled circle 
marks where the final relaxation time r = lO^'^/F is read off, <I>(t = r) = 0.1). Curve H lies at 
the critical point e = for 7/F = 10"^'^, while I to M keep 7/F = 10^^'^ but increase e = 4"ei, 
n = 0, . . . , 4. The horizontal bars indicate the heights of the non-ergodicity parameters / for e > 
(compare curves H to M), which would be approached at long times for 7 = 0, ^{t — > 00) = /. 

A. Transient fluctuations 

Figure |^ shows correlators of the F^'^^ model for various shear-rates and distances from 
the non-ergodicity transition at e = 0. Curve G there corresponds to a non-sheared fluid 
state close to the transition, e = —ei = — 10~'^''^^ and 7 = 0. This shows the typical two-step 
relaxation pattern with microscopic short-time dynamics for tT = 0{1), followed by the 
approach to an intermediate plateau at fc ~ 0.293, and the final decay characterised by the 
(final or a-) relaxation time r. For positive separation parameters, e > 0, the final decay is 
absent and the correlators approach finite long time limits $(t 00) = f > fc (not shown; 
these /'s are indicated by horizontal bars at the left of the figure). Including a finite shear 
rate corresponding to a small but finite bare Peclet number Peo = 7/F in the model, little 
effect of shear on the short time dynamics is seen because of Peo <^ 1. A drastic effect on the 
final decay, however, is seen in curves A to E, because the dressed Peclet (or Weissenberg) 
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number Pe= 7r is not negligible. Moreover, all glassy curves {e > 0), which would stay 
arrested for 7 = 0, are seen to decay by a shear-induced process whose time scale is set by 
the inverse shear rate, and whose amplitude depends on the distance to the transition. 

The stability analysis of Eqs. (0,^ describes the correlators for a window around $q ^ 
and a finite window in time, which both can be estimated from the condition \G{t)\ ^ 1; for 
7 = these have been worked out in detailii. As the analysis in Appendix B shows, (7^)^ 
dominates for long times in Eq. (|^), and therefore always 

^(t - 00) ^ -t/rW , with r(^) = \/^^^. (12) 

Hence, under the imposed shear flow, density fluctuations always decay, as this decrease 
of Q{t) for long times initiates the final relaxation of $q(t) to zero. (In this region the 
corrections to Eq. (^ of higher order in Q become important.) Arbitrarily small steady 
shear rates 7 melt the glass, as expected, and for small Pcq there appears a separation of 
time scales between short-time motion and the shear induced final decay; see curves H to M 
in figure ||. 

For more detailed insight into the shapes of the relaxation curves, we turn to the ISHSM. 
Figure ^ displays density correlators at two densities, just below (panel (a)) and just above 
(panel (b)) the transition, for varying shear rates. Again, in almost all cases Pcq is negligibly 
small and the short-time dynamics is not affected. In the fluid case, the final or a-relaxation 
is also not affected for the two smallest Pe values, but for larger Pe the non-exponential a- 
relaxation becomes faster and less stretched; see the inset of fig. |5(a). (While the techniques 
of Appendix B allow us to discuss this shape change, we do not do so here for lack of space.) 
The glassy curves at e > 0, panel (b), exhibit a shift of the final relaxation with r^"^) 
and asymptotically approach a scaling function $+(t/r'^^)). The master equation for the 
"yielding" scaling functions in the ISHSM can be obtained from eliminating the short- 
time dynamics in Eq. (^. After a partial integration, the equation with ^q{t) = is solved 
by the scaling functions: 

= <W - i J^di' m+(t- ?)$,+ (?) , (13) 
where i = t/T^'^\ and the memory kernel is given by 

^tii) = ^ E (i) '^tii) ■ (14) 

k 
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FIG. 3: Normalized transient density correlators ^q{t) of the ISHSM at wavevector q = 3.4/d below 
(panel (a) at (?!) = — 10~^) and above (panel (b) at ^ = + 10~^) the transition for increasing 
shear rates Pcq = 9"^ * 10~^ with n = 0, . . . , 10 from right to left; the distances correspond to 
£ = ±10~^'^^. Curves for ra = 9, 10 carry short and for n = 8 long dashes; note the collapse 
of the two short dashed curves in (a). The insets show the data rescaled so as to coincide at 
$(t = r) = 0.1. 
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While the vertex is evaluated at fixed shear rate, 7 = y I^ItT' depends on the equilibrium 
parameters. The corresponding results for the Fi^^ can easily be obtained. 

The importance of the /5-correlator Q for the derivation of the above scaling functions is 
because Q, according to Eq. (0), depends on r^^^ and not on internal time scales like r. 
It provides the initial conditions for Eq. ([13|) , which follow from the analysis of Eq. in 
Appendix B. In the glass, e > 0, Eq. (pSf) leads to 



<^Ui^o) = r + h, 



i-——p+o{-.m 



(15) 



and for somewhat longer times, i ^ y/e, according to Eq. (p9D, this merges into: 

<l>+(v^ « i « 1) = ~hj. (16) 

At the transition, e = 0, Eq. ( P8| ) shows that Eq. ([T6|) holds down to t ^ 0. The 
scaling functions defined thereby describe the yielding behavior of a sheared glass where the 
relaxation time is set by the external shear rate and internal relaxation mechanisms are not 
active. 

Figure ^ shows the transient shear modulus (7(^,7) of the ISHSM which determines the 
viscosity via 77 = J^dt g{t,'y). It is the time derivative of the shear stress growth function 
r]^{t,'j) (or transient start up viscosity; here, the labels the shear history^), g{t,'j) = 
^r7+(t, 7), and in the Newtonian- regime reduces to the time dependent shear modulus, 
G{t) = g(t,'y = 0). The git,^) shows all the features exhibited by the correlator of the 
schematic model, and by the density correlators in the ISHSM, and thus the discussion based 
upon Q and the yielding scaling law carries over to it. But in contrast to the correlator of 
the Fij\ and more so than the <l>g, the function g{t,'j) becomes negative (oscillatory) in the 
final approach towards zero, an effect more marked at high Pe. This behavior originates in 
the general expression for g{t,'y), Eq. (|1|), where the vertex reduces to a positive function 
(complete square) only in the absence of shear advection. A overshoot and oscillatory 
approach of the start up viscosity to the steady state value, ri^{t — * 00, 7) — » rj^^y), therefore 
are generic features predicted from our approach. 

In the discussion that follows, the leading corrections to the scaling law for the yielding 
correlators will be required. The analysis of Eq. (^ suggests that to next order for £ > 0: 

%{t) ^ $+(t) + 6{e, |7to|) $+(t) + 0{6') , (17) 
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FIG. 4: Transient non-Newtonian shear modulus g{t,^) of the ISHSM in units of kBT/d? for the 
packing fractions (j) = (pc^ 10~^ (e = ±10"^'^^; dashed/ dotted hues, respectively) for increasing 
shear rates Peo = 9" * 10~^ with n = 0, . . . , 8 from right to left; note the collapse of fluid lines for 
the smallest Peo- The solid line gives 5(^,7) for (j) = (pc and Peg = 10~^. The inset shows the data 
rescaled so as to coincide at g{t = r, 7) = 5; note the collapse of the n = 6,7 & 8 curves for both 
(j) > (pc and (j) < (pc- 

where to is a matching time, and the correction exhibits a weak divergence for short 

rescaled times, ^ 0) oc t-(2A-i)/A^ follows from Eqs. (|B|,^9D. The inset of fig. ||(b) 



shows the rise of the correlators above the yielding master function at short times and for 
very small Pcq. Except for the fact that it is integrable, this function is of no further interest 
here. The small parameter 6, however, sets the scale of the corrections and (as shown in 
Eqs. ( P8| , P^ )) exhibits the following scaling properties 

f |7tor with m = for e = 

5 oc r ' ' ^ (18) 

[ e™' with m' = for |7to|^ < £ < 1 

These results will be used below. 



B. Flow curves 

In our approach, steady state properties of the sheared system are obtained via time 
integrals — from switching on the rheometer at t = up to very late times when the 
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system has relaxed into the non-equihbrium stationary state. The evolution of the system 
is approximated by following the transient density fluctuations, which were discussed in the 
previous section for two specific simplified models. Equation (|1]) gives the thermodynamic 
shear stress in our fully microscopic approach, while Eqs. (0,|TD|) give simplified expressions 
for it using the two models. Inserting the transient density correlators from section |1V A| , 
the stress versus strain rate curves ( "flow curves" ) of the two models can be discussed. Such 
relations o'{'y) often are postulated en route to deriving phenomenological "constitutive 
equations" for nonlinear flow behaviour, whereas our approach leads, in principle, to their 
microscopic derivation. Figures H to || present the results, the first two as plots of a versus 
7, while the second pair shows r] versus 7. We will discuss the general qualitative features 
jointly for both models. 




FIG. 5: Steady state "stress" a = 7(r) in the schematic F^2^ model as function of Peo = 7/r at 
various distances from the transition: e = (solid line) and e = ib4"ei for n = —1,0, ..,4 (while 
El = lO^'^'''^) with short (long) dashes for fluid e < (glassy e > 0) curves. At e = 0, a dotted 
line gives a = cr+(l + 2.48 7'''^^^) for 7^0, with 0"+ = 0.10, while a thin dot-dashed line gives 
(Tc = 7 for 7 — > 00. The inset shows the finite limits cr+ = a{e > 0, 7 ^ 0+) obtained above the 
transition; the dot-dashed curve is = + 1.88-^6. 

In the fluid, we find a linear or Newtonian regime in the limit 7^0, where we recover 
the standard MCT approximations for the transport coefficients of a viscoelastic fluid§§. 
Hence a 7^7(7 = 0) holds for Pe -C 1, as shown in figure where Pe calculated with 
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the structural relaxation time r is included. As discussed in section |T|, the growth of r 
(asymptotically) dominates all transport coefficients of the colloidal suspension and causes 
an proportional increase in the viscosity 77. For Pe > 1, the non-linear viscosity shear thins, 
and cr increases sublinearly with 7. Without analysing these complicated flow curves in 
detail here, we note that additional, shorter time scales than r enter; these cause the shape 
change of the density correlators shown in figure ^(a) and also affect a and rj. The stress 
versus strain rate figures |^, ^ clearly exhibit a broad crossover between the linear Newtonian 
and a much weaker (asymptotically) 7-independent variation of the stress. 

Replotting the identical data as viscosity versus strain rate (figures 0, ||) these subtle 
features get compressed by the requirement of a much larger scale range on the viscosity 
axis; hence plotting stresses should prove more telling, close to a glass transition, than 
plotting viscosities. The latter can reveal subtle non-asymptotic corrections, but only if 
they are replotted on a smaller scale as done in the insets of figures |^, ^. There, various 
apparent power-law exponents could be fitted to the numerical data. Considering the ISHSM 
as a model for colloidal suspensions, the high-shear limiting viscosity contribution needs to 
be included; neglecting hydrodynamic interactions, this is set by the solvent viscosity rjoo- 
The corresponding dot-dashed curves in figures ^ |] show that a rather close approach to 
the critical packing fraction is required for the structural stress captured in Eq. (|^) to 
dominate. Considering that hydrodynamic interactions cause an appreciable increase of the 
high-shear limiting viscosity over the solvent onell, this condition appears severe. However, 
it is obviously satisfied for systems that are nonergodic at rest (e > 0), if small enough 7 
can be achieved. 

Above the transition, the quiescent system forms an (idealized) glassEl'0 which exhibits 
finite elastic constants. The transverse elastic constant Goo describes the (zero-frequency) 
Hookian response of the amorphous solid to a small applied shear strain 7, so that a = Goo7 
for 7 — > 0. If steady flow is imposed on the system, however, we find that the glass yields 
and is shear melted by arbitrarily small shear rates. This fiuidization is not simply a trivial 
consequence of advection of the particles with the flow (Taylor dispersion is included in 
neither model) but implies that particles are freed from their cages and diffusion perpen- 
dicular to the shear plane also becomes possible. Hence, within our approach, infinitesimal 
steady shear leads to true melting of the glass and not merely plastic flow of it. Any finite 
shear rate, however small, sets a finite longest relaxation time, beyond which ergodicity is 
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FIG. 6: Steady state shear stress a in units of ksT/d^ versus Peo = '^(f'/D^, for the ISHSM at 
various distances from its glass transition, (/) — (/)c as labeled; circles correspond to fluid, diamonds 
to the critical, and squares to glassy densities. For the fluid cases, (j) < (f)c, dashed lines indicate 
Newtonian fluid behavior, a = rij, while vertical bars mark Pe= 7r = 1, with the structural 
relaxation time taken from ^g^j/dit = t) = 0.1. The stress which would additionally arise from 
the background solvent viscosity, a = 7r/oo) is marked by a dot-dashed line. For the critical 
density, (pc, the critical yield stress, = 6.04, is shown by a horizontal bar, and the dotted line 
a = a'^{l + 0.897°'^^^) holds for 7^0. Crosses connected by a thick dashed line give a at (p = (pc 
but for a different numerical discretization of the memory kernel (3-times larger cut-off /cmax)- The 
inset shows the rise of the dynamical yield stress = a{e > 0,7 ^ 0-|-) in the glass together with 
the fitted square root asymptote, = + 112^/(j) — cpf.- 

restored; see the discussion of figures 0, ^. 

Nonetheless, a finite limiting stress (yield stress) must be overcome in order to maintain 
the flow of the glass: cr(7,£: > 0) ^ ^"^(e) for 7^0. To understand this better, note that 
for £ > and 7 — > 0, the time r*^'^-' for the flnal decay , Eq. (|T2|), can become arbitrarily 
slow compared to the time characterizing the decay onto fq. Inserting the scaling functions 
from Eqs. ([TsD to (p^Bf ) into the expressions Eqs. (|7|,[lOD for the stress, the long time 
contributions separate out. Importantly, the integrands containing the functions depend 
on time only via 7^, so that nontrivial limits for the stationary stress follow in the limit 
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7^0. In the ISHSM for e >0, this is given by (for 7 > 0): 



a 



607r2 



fdijdkk'^^ , (19) 

k{t) 

where k{i) = k (^1 + (7^)^/8^ , and the fixed reduced shear rate 7 was defined after Eq. 
(P^). The corresponding result in the ^^^2'* model is simply = 7 J^di $^(t). 

The existence of a dynamic yield stress in the glass phase is thus seen to arise from the 
scaling law in Eqs. (|l3| , p^ , whose decay is initiated by the shear-induced asymptote of the 
correlator Q in Eq. ([121) . Both the ISHSM and the schematic model share this feature. The 
yield stress arises from those fluctuations which require the presence of shearing to avoid their 
arrest. Importantly, even though requires the solution of dynamical equations, within 
our approach it is solely and uniquely determined by the equilibrium structure factor Sg at 
the transition point. Assuming the connections of MCT to the potential energy paradigm 
for glasses, as recently discussedHS, one might argue that cr+ arises because the external 
driving allows the system to overcome energy barriers so that different metastable states 
can be reached. This interpretation would agree with the recent spin-glasJ^^ and soft-glassy 
rheology0iB§i approaches. Our microscopic approach indicates how shear achieves this in 
the case of colloidal suspensions. It pushes fluctuations to shorter wavelengths where smaller 
particle rearrangements cause their decay. 



The increase of the amplitude of the yielding master functions see Eq. (15) originates 
in the increase of the arrested structure in the unsheared glass (via the non-ergodicity 
parameters fq — fg oc y/e). In consequence the yield stress should rapidly increase as one 
moves further into the glass phase — oc y/e should (approximately) hold. Indeed, 
the insets of figures |, ^ show a good fit of this square root increase to the numerical data. 
Our ability to make such a prediction is highly significant: the yield stress is an important 
non-linear property of the arrested state itself, it can be calculated here, because the yield 
stress is governed by the onset of melting under shear, which is itself a glass transition — 
not at (f) (pc but at 7 — > 0+. 

An intriguing power-law increase of the stress above the yield value at the critical point 
was noticed in Ref.lii and is also included (as dotted lines) in figures ^ and |^. It results from 
the leading corrections to the yielding scaling law summarized in Eqs. (|T^|IB]). Inserting 
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FIG. 7: Shear dependent "viscosity" rj = (r) (in units of 1/r) in the schematic F^'^'^ model as 
function of Pcq = 7/r at various distances from the transition (data and hne styles as in figure |^ . 
The inset shows the fluid curves e < rescaled with the viscosity r/o = ^(7 = 0) of the unsheared 
situation. Two power-laws (dashed lines) indicate the minimal and maximal slopes of the curves 
around 7770 = 10^. 

those expressions into the integrals for the stress leads to the small shear rate expansion 

a = (T+ + |7tor CI + 0(17^01"") , for |£| < |7to|^ , (20) 



where the constant a is given by an identical integral to that in Eq. ( [T9| ) but with 
replaced with the correction <!>"'". For the ISHSM the exponent m turns out as m = 0.152 
while it is m = 0.143 for the model. As this prediction is based on the universal stability 
equation of a yielding glass, Eqs. (0,|]), and on the existence of the yielding scaling law, Eqs. 
(|I^,0), our approach suggest that this Hershel-Bulkeley0 flow curve may hold universally 
at a glass transition point, with the exponent m depending on the material via the static 
structure factor Sq. 

Another prediction from our asymptotic analysis of the leading corrections is borne out 
in the numerical calculations. Deep in the glass, Eqs. ([T7| , [T8|) and (|20|) state that the leading 
correction to a is independent on the shear rate 7: a = + e^' + 0{e^'^'\ This explains 
why in figures ^ and |^ the stress starts to rise with 7 appreciably only when Fcq becomes 
non-negligible. 
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FIG. 8: Steady state non-Newtonian shear viscosity r] in units of ksT / D^d versus Peg = ^(f / Dq, 
for the ISHSM at various distances from its glass transition, (p — (pc as labelled; data and line styles 
from figure |6| (the dot-dashed horizontal line gives the neglected background solvent viscosity). 
The inset shows the fluid curves e < rescaled with the viscosity rjQ = 77(7 = 0) of the quiescent 
fluid. Two power-laws (dashed lines) indicate the minimal and maximal slopes of the curves around 
7% = lO^. 

For the ISHSM, a caveat concerning the numerical results is required for shear rates 
beyond about Pcq > 10^^. The wavevector integrals in Eq. (|^) do not converge properly for 
large k, and the results thus depend on the chosen cutoff, fcmax = 39.8/d. This is implied 
by the crosses included in figure which were calculated for (p = (p^ but larger cutoff, 
kmax = 119.8/(i. While only small differences (between the crosses and diamonds) remain for 
Peo < 10~^, for larger shear rates the more accurate integrals lead to larger stresses. While 
this could be an artefact of the ISHSM, figure § shows that the time-dependent non-linear 
shear modulus g{t, 7) decays rapidly for such large Peclet numbers. Because its initial value, 
g{t = 0, 7), is the instantaneous modulus of the unsheared system, g(t = 0, 7) = G(t = 0) = 
{o'ly)/V, this depends strongly on the particle interaction potential, and, for hard spheres 
should actually divergeiiS. As changing fcmax strongly influences G{t = 0), we presume 
that finite /cmax effectively softens the particle repulsion. Only calculations for more general 
potentials can show whether this indicates a strong dependence of the non-linear stress on the 
steepness of the potential for not very small Peclet numbers. Interestingly, the glassy yield 
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stress does not exhibit this strong dependence. It should however be noted that the ISHSM 
underestimates the effects of shearing as the ratio cr+/G^ = 0.33, is overestimated0ll0. 

C. Structural versus non-structural features 

Recently, computer simulations of sheared atomic glass formerJ^ljiiii have been per- 
formed. While our microscopic theory is firmly based on assuming colloidal dynamics, in 
the simplified models, the effects of a different short time dynamics can easily be studied. 
Without shear, the MCT has found that the long-time structural relaxation is independent 
of the microscopic short time dynamicsBil. (The latter only sets the overall time scale via 
the matching time to-) In the F^2'' model, in order to mimic Newtonian dynamics, Eq. (^ 
can be replaced by 

$(t) + u <i>(t) + |$(t) + Jdt' m{t - t') <i>(t')| = , (21) 

while Eq. (^ remains. Here, is a microscopic vibrational frequency and v a bare damping 
coefficient. Varying v/Vl shifts to? and the thin lines in figure |^ show the effect on the stress 
of the critical glassy state. The factor x = {t^^Q) /{tQ^T) varies between 0.3 and 5. As 
discussed in the context of the yielding scaling law, Eqs. (0,0) and Eq. (0), the yield 



stress o""*" is, according to our approach, a purely structural property which is independent of 
the microscopic transient dynamics. Beyond the limit 7 — > 0, the one microscopic matching 
time to enters the structural dynamics and affects the rise of a. The data can be collapsed 
over a substantial window by replotting them versus jto as explained by Eq. (^OD; see the 
bold lines in figure ^ This suggests that the expansion Eq. ( pO| ) could be extended to 
higher orders and might then obtain a larger range of validity. As the average relaxation 
time To = (t(7 = 0)) of an unsheared fluid state close to the transition is proportional to to, 
replotting the rescaled relaxation times ("viscosities"), (r(7))/ro, versus rescaled shear rate, 
7x0, ehminates the to dependence and leads to data collapse even though the the expected 
asymptotic power law, 77 oc 7"^ is still strongly distorted. 



To conclude, as argued already in section |IV B|, even under shearing the structural relax- 



ation remains independent of short time microscopic effects, except for shifts in the matching 
time to- According to this, simulations with either Newtonian or Brownian dynamics (that 
are otherwise identical) should thus observe the same nonlinear rheology near the glass 
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transition. Rescaling with to eliminates transient features and the nonhnear stress is only 
determined by the static structure factor for Peo ^ 1. This last statement bears an im- 
portant caveat though. The time is rather fast and thus not well separated from other 
microscopic time scales, so that corrections from finite Peg need to be anticipated. This 
remains a task for future work. 




FIG. 9: In the main panel, steady state "stress", a = 7(t), versus rescaled Peo for the schematic 

(7) 

model with different short time dynamics at the critical point e = 0; a dotted curve gives the 
Hershel-Bulkeley law from figure ^. In the inset, shear rate dependent "viscosity", rj = (r), for 
different short time dynamics below the transition, e = —4?ei = — 10~^'^^, rescaled with the shear 
rate independent 770 = 77(7 = 0); a dotted lines gives a power-law fit. The viscosities rjo at this e 
also determine the rescaling factors used in the main panel: /t^^ = Vo^ /Vo^- Solid lines give 
Brownian dynamics results, while dashed curves correspond to Newtonian dynamics with z^/il = 5, 
1 (both short), 0.5, 0.1 and 0.01 (all three long dashes). Thin lines in the main panel give the 
unshifted a for = 5 and 0.01. 

D. Non-linear Maxwell model 

In 1863, Maxwell suggested a simple model for the linear rheological properties of a 
viscoelastic fluid, which has been the cornerstone of the phenomenology of glassy systems 
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since then. He suggested that the time dependent shear modulus decays exponentially, 
G{t) = Goo exp —t/r, where Goo is the solid shear modulus and r the structural relaxation 
time. Viscosity and stress follow as usual, cr = 777 = 7 J^dt G(t). Our results suggest, as 
a toy model, a non-linear extension of Maxwell's model which captures the competition of 
structural arrest and shear-induced motion. It is obtained simply by postulating that the 
transient shear modulus can relax via two independent processes, an internal one charac- 
terized by r, and an induced or driven one characterized by r^^^ = c*/|7| (with some 
numerical constant): 

^^-"(t,7)=G.oe-/^^ with ^ = 1 + :^. (22) 

From this the expression for the flow curve follows immediately as 

(23) 

which exhibits a simplified scenario of the yielding-to flow transition which we derived within 
our microscopic approach; a glass would correspond to 1/r = 0, in which case the flow curve 
is that of a Bingham plastic. That toy model is not totally trivial, in that it gives a non- 
analytic dependence on the shear rate. 

V. CONCLUSIONS AND OUTLOOK 

Building on Ref.El, we have presented a microscopic theory of the nonlinear rheology of 
colloidal fluids and glasses under steady shear. This predicts a universal transition between 
shear-thinning fluid flow, with diverging viscosity upon increasing the interactions, and solid 
yielding, with a yield stress that is finite at (and beyond) the glass point. Numerical calcu- 
lations presented here within progressively more simplified models support these universal 
predictions. The novel yielding behaviour is seen to arise from a competition between a 
collective caging of particles driven by increasing local order, and shear advection of the in- 
volved density fluctuations to smaller wavelengths, where Brownian particle motion relaxes 
them more effectively. 

The shear-melting scenario can be rationalized solely from the stability analysis of a 
yielding solid summarized by Eqs. (||,|D. As we tried to work out carefully, the properties 
of the correlator Q and some (natural but as yet unproven, except for the models discussed 



(7 = 7 



Voo + 



GnoT 



1 + |7|r/c. 
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here) assumptions about its extension to longer times suffice to explain our major predictions. 
These include: (i) The existence of a dynamic yield stress of glasses, which increases strongly 
with increasing density beyond the transition, (ii) The relatively small increase of the 
potential part of the stress in the solid above the yield value upon increasing 7. {iii) The 
prediction of a Hershel-Bulkeley flow curve at the critical point, with a power law index 
fixed by the quiescent static structure factor, (iv) An asymptotic power-law shear thinning 
viscosity, rj oc 7"^ in the fluid (this follows from the existence of a yield stress), {v) The 
observation that motions for times shorter than the final relaxation time r influence the 
stress or viscosity beyond Pe = 7r <^ 1. This finding, discussed for a schematic model in 
section IV C , predicts strong sub-leading corrections to points {n)-{iv) above and shows that 
more detailed calculations (which ought to include hydrodynamic interactions) are required, 
especially when no clear separation of time scales is present. 

Reassuringly, the central non-linear stability equation, Eq. (H), which lies at the core 
of our predictions, can be argued on more general grounds than the explicit derivation 
we gave. Without shear, its predictions have been tested experimentallylliSill, and its 
structure thus appears established (with some confidence certainly for e < 0). The simplest 
inclusion of shear naturally leads to a term — (7t)^ as the essential perturbation for long 
times. (Obviously the sign of 7 must not matter, and the term must be negative to induce 
relaxation.) Also, because Eq. (|^) without shear is a scaling equation without inherent 
time-scale (viz. ^(y * t, 7 = 0) is a solution for any ?/ > if it is for y = 1), the shear- 
rate 7 can only enter multiplied by time t itself. Thus we expect Eq. (^ to hold within a 
wider context of non-ergodicity transitions and hope that the predictions {i) to (v) above 
may apply to a broader class of scenarios describing yielding of metastable solids and shear 
thinning of viscoelastic fluids. 

On the other hand, the approach we have outlined cannot address the shear thicken- 
ing behaviour that, for many colloidal materials, occurs at higher flow rates than those 
addressed here (but in some cases appears to preclude shear thinning altogether). It has 
been suggestedS^ that this jamming phenomenon is connected with a stress-induced glass 
transition. This avenue will be explored in future work on a version of the schematic model 
in which Eq. (^ is modified to include explicit stress- (as well as strain-rate-) dependence^. 
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APPENDIX A: BIFURCATION ANALYSIS OF THE MICROSCOPIC AP- 
PROACH 

To determine steady state quantities hke the stress, the transient density fluctuations 
must be found within the suggested microscopic approximation scheme. With t denoting 
the time after switching on the shear field, the normalized density correlators obey the exact 
equation of motion: 

+ T^{t) |<l>q(t) + Jdt' m^it.t') •i'qlt')} = , (Al) 

where the "initial decay rate" becomes time-dependent with shear and exhibits Taylor disper- 
sion, T^{t) Sq = Do {(f + qxQyit + {qxqyit + ql'^H'^)Sq - {qxqy/q)iS'^. The memory func- 
tion mf^{t, t') can only be found after approximations, and, in the MOT spirit of our approach, 
is approximated by projecting the fluctuating forces onto density pairs and factorizing the 
resulting pair-density correlation functions: 

^ ^EK;S(i,0 Mt-t')%-u{t-t') . (A2) 

k 

The vertex V"^'*'-' is evaluated in the limit Peo <^ 1 but for large times so that jt and jt' are 
finite. As 7 0, it reduces to the standard MOT vertexlll and like the latter is uniquely 
determined by the equilibrium structure factor, Sq. For long times, it vanishes. 

The (approximate) microscopic equations contain the central bifurcation singularity of 
MCT whose universal properties can be obtained from a stability analysis close to the 
singularity. It is these universal properties which the simplified models share. Around a 
quasi-arrested glassy structure, characterized by parameters f^, the Eqs. ( |A1| , [A2D can be 
solved with the ansatz 

$qW = /.^+(l-/,T^qW, (A3) 

where we use the fact that without shear the system is isotropic and the (critical) glass 
form factor therefore depends on the magnitude q only. Splitting the vertex close to the 
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transition into its critical value and deviations 



(t, t') = y + 6v:s. + (t^ t') m 



the stability of can be determined. Introducing the (standard) abbreviations for time 
independent values, = 2^ X^k ^^^^fkfp ^'^d £q = jj^ Z^k ^Ki,k/fc/p) ^-^id for matrix coeffi- 
cients in the expansion in G, = ^ Ek' Ki!k'Vfc(l - fpfh,k' and = Ek' 'KJ^'^'II - 



- /p)^5fc,fc'5p,|q-k'|, Eq. (^) for long times becomes 



fc 

+ 



£,-^EkAKi,^(t,0)/,%=- 

Efcp + I Efc - /orft' ^.(t - t')g,{t') + ... . (A5) 

Here, a joint factor [1 — f^) was divided out and (. . .) denotes both terms of higher order in Q 
and those connected to the short time dynamics which become negligible for long times. For 
large enough vertices, the first two terms in Eq. (|A5| ) of order have a finite solution 
which defines (idealized) glass points. A critical (glass) point, viz. bifurcation singularity in 
Eqs. ( |70]JA2|) , lies where the square bracket in Eq. ( |A5|) vanishes because has eigenvalue 
1, and the linearized equation cannot be solved for Q{t), which gives the dynamics along 
the unstable direction, Qq{t) = Q{t)eq. Here the eigenvectors of are denoted as and e^, 
which satisfy X^fc Cqk^k = and Y.q ^qC^^ = ^fc and the (conveniently chosen) normalizations 
Y,q (iq^q = Hq ^qi)- ~ fq)^q^q = 1- The Solvability equation for Eq. (|A5|) results from requiring 
the contributions in the second and third lines to be perpendicular to iq, and without shear 
Eq. d^) follows upon the definitions: e = J2q€-q£q and A = J^qkp^qCqkp^k^p; for more details 
se With shear, a time dependent contribution arises from AVq^J(t,0), and because 

of spatial isotropy in the limit of small shear, this reduces to: 

^^'^ = li^fj, ^ E ^q AKi^J (t, 0) /|Vk| , ( A6) 

where Cq = ^J^2y takes care of three-dimensional integration. For colloidal hard spheres 
described with the Percus-Yevick structure factoJ^, Monte Carlo integration estimates give 
c(^) ^ 3. 
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APPENDIX B: ANALYSIS OF SCALING EQUATION (3) 



The so-called /3-correlator Q is the solution of Eq. with the prescribed short time 
behavior Q{t ^ 0) = {to/t)^ where to is a matching time to the short time transient (deter- 
mined by r), and the "critical" exponent a obeys A = r^(l — a) /r(l — 2a). The derivation of 
this behaviour, which solves Eq. (^) at the singularity for long times, and the properties of 



Q without shear can be found in Refs.ES^^. With shear, the analysis of Eq. (|^) resembles 
the one of the extended MCTil, where a relaxation kernel describing thermally activated 
hopping motion lead to a term linear in t instead of the quadratic 7^t^. Because, in the main 
text, we focus on the yielding behavior of glassy solids, a discussion of the consequences of 
this quadratic term close to the glass transition and in the glass shall now be given for the 
regime where shear-induced motion dominates. 

A useful aspect of Q are its homogeneity properties. With i7 > an arbitrary scale, it 
obeys 

g{t, £, 7) = ^ [i = Q{t/to),e = en-^\ ^ = 7^0^^"^'+"^) , (Bl) 

and thus does not change shape if the two control parameters are varied on a scaling line 
7^0 = 'y{e/e)''^^"'^^^'^"'\ Moreover, this property enables one to define three regimes, (i) 
e <^ —Ea/ far below the transition in the fluid region, where finite shear barely distorts the 
fluid /^-correlator; (ii) \e\ <^ e^, in the transition region, where the dynamical anomalies 
are most pronounced; and [iii) e ^ far above the transition in the shear yielding glassy 
state. Here, e-y = l^ytoP"^''^^'^^ defines a natural scale for the separation parameter e, and 
throughout the following, <^ 1 and 7^0 -C 1 is assumed. 

Equation (Rf) can then be solved by power series expansions with the general ansatz 



g = t^ J2 ^nt'''' , (B2) 

n=0 



which yields 



_ ^(7) 



£ - (7 tf = a^ro,ot'" + 2aoairo,it'"+'' + E E r„_„>.a„_„,a„. . (B3) 

n=2 n'=0 



Here, an abbreviation employing Gamma functions r(a;) is used: 



_ T{l+u + nv)T{l + u + n'v) _ 

r(l + 2M + (n + nO^) ^ ' 
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1. Glass for intermediate times 



For any 7, a solution to Eq. (|B3|) can be found for e > by choosing u = 0. This requires 
ao = \J^/{.^ — A), V = 2 and ai = —0*^^)7^/(20:0(1 — A)), while the higher coefficients can be 
found from a simple recursion relation. The result (with = e^/^^"^)) takes the form 

c(^) , t 



^(e > 0, tf < t < 4) = ^ ( ^e^, +1, ^toe-^ ) = Git/t 



2 

(B5) 

where = a/^/|7| is a time scale that follows from balancing the first two terms in Eq. 



(0). Without shear, the correlator would adhere to ye/{l — X) for long times in the glass, 
while for finite 7 it relaxes below this value. The series Eq. (p5| ) gives the solution of Eq. 
(0) for a finite window in time only, tf ^ t ^ ti,, as can be seen from its failure to match 
to the {to/t)"- at short times, and as it has a finite radius of convergence of the order of tf,; 
the latter property can easily be determined from the recursion relation for the an>2- The 
obtained series thus describes, to order ^/e, the initial stage of the final decay down from 
f, = r, + ^e/{l~X). 



2. Transition region for long times 



If, close to the transition, e = can be set (region (ii)), then the choice u = 1 requires 
Oo = \Jc^'^V{.^ " \) \i\ in Eq. (p3| ). To match the resulting long-time series, higher a„ are 
needed, and can be recursively be found, iff Vq i vanishes: 

r(2)r(2 + t;) 



To,: 



A = 



^ V 



(B6) 



T{3 + v) 

Obviously, coefficient ai then cannot be determined from Eq. (P3|), but needs to be found 
from matching to Q at shorter times. Its scaling with 7 can easily be obtained from the 
homogeneity statement, Eq. ( pll) : ai = aoai^^^^ , where ai is a matching constant (pre- 
sumably of order unity), and the shear-rate dependent /5-scaling time appears 



t-y = to |7^o| ■ 

As announced previously, the scaling function Q then follows (with f2 = |7to|~ 

i7toi^ G{tm) -- 



(B7) 



£-g -,0,1 



7(7) 



-\i\t 



ai — ^ + 



(B8) 
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It matches to the decay of from unity down to f^, which determines ai, and describes 
the initial stage of the shear-induced decay of the correlators down from to zero. 



3. Glass for long times 



Deep in the glass, region [iii), the same reasoning as in the case just studied can be used 
to determine ai from matching — now to the series in Eq. (|B5D. With u = 1 and v from Eq. 
([B6D, the scaling of the undetermined coefficient is fixed by requiring it to extend Eq. 



from around tf, to longer times: ai = ao&'itfj^^^ , where &[ is another matching constant. 
The resulting series continues the shear induced decay of the correlators from the plateau 
value, as follows: 

^(^) . . r ,t 



\ito\^+^ ^ (^,— ,1 



^ Q{tlh) 



-\i\t 



, , - , 2A-1 



(B9) 
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